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The stochastic Schrodinger equation, of classical or quantum type, allows to describe open 
quantum systems under measurement in continuous time. In this paper we review the link between 
these two descriptions and we study the properties of the output of the measurement. For simpli- 
city we deal only with the diffusive case. Firstly, we discuss the quantum stochastic Schrodinger 
equation, which is based on quantum stochastic calculus, and we show how to transform it into the 
classical stochastic Schrodinger equation by diagonalization of suitable quantum observables, based 
on the isomorphism between Fock space and Wiener space. Then, we give the a posteriori state, 
the conditional system state at time t given the output up to that time and we link its evolution to 
the classical stochastic Schrodinger equation. Finally, we study the output of the continuous meas- 
urement, which is a stochastic process with probability distribution given by the rules of quantum 
mechanics. When the output process is stationary, at least in the long run, the spectrum of the 
process can be introduced and its properties studied. In particular we show how the Heisenberg 
uncertainty relations give rise to characteristic bounds on the possible spectra and we discuss how 
this is related to the typical quantum phenomenon of squeezing. We use a simple quantum system, 
a two-level atom stimulated by a laser, to discuss the differences between homodyne and hetero- 
dyne detection and to explicitly show squeezing and anti-squeezing and the Mollow triplet in the 
fluorescence spectrum. 



I. INTRODUCTION 

A big achievement in the 70's-80's was to show that, 
inside the axiomatic formulation of quantum mechan- 
ics, based on positive operator valued measures and in- 
struments [H E] , a consistent formulation of the theory 
of measurements in continuous time {quantum continu- 
ous measurements) was possible [2HS]- Starting from 
the 80's, two other very flexible and powerful formula- 
tions of continuous measurement theory were developed. 
The first one is often referred as quantum trajectory the- 
ory and it is based on the the stochastic Schrodinger 
equation (SSE), a stochastic differential equation of clas- 
sical type (commuting noises, Ito calculus) p, [71 [5HT7]. 
The second formulation is based on quantum stochastic 
calculus |18H20j and the quantum SSE (non commut- 
ing noises, Bose fields, Hudson-Parthasarathy equation) 
[Ml li [II mi nil nil Hi]. The main applications of 
quantum continuous measurements are in the photon de- 
tection theory in quantum optics {direct, heterodyne, ho- 
modyne detection) [QHTtI [2Tti25] . While the classical SSE 
gives a differential description of the joint evolution of the 
observed signal and of the measured system, in agreement 
with the axiomatic formulation of quantum mechanics, 
the quantum SSE gives a dilation of the measurement 
process, explicitly introducing an environment which in- 
teracts with the system and mediates the observations. 

In this paper we start by giving a short presenta- 
tion of continuous measurement theory based on the 



* Also at Istituto Nazionalc di Fisica Nucleare, Sezione di Milano, 
and Istituto Nazionale di Alta Matematica, GNAMPA 



quantum SSE (Sees. |TT] and III). We consider only the 
type of observables relevant for the description of homo- 
dyne/heterodyne detection and we make the mathemat- 
ical simplification of introducing only bounded operators 
on the Hilbert space of the quantum system of interest 
and a finite number of noises; for the case of unbounded 
operators see [33H1S]. 

In Sec, mil we show how to derive the classical SSE and 
the related stochastic master equation (SME). The key 
point in the step from the quantum SSE to the classical 
SSE is the introduction of an Hilbert space isomorphism 
which diagonalizes a suitable complete set of quantum 
observables. The classical SSE and the SME give both 
the probability distribution for the observed output and 
the a posteriori state, the conditional system state given 
a realization of the output. These equations are driven 
by classical noises, but, in spite of this, they are fully 
quantum as they are equivalent to the formulation of con- 
tinuous measurements based on quantum fields. It is just 
in this formulation that the probabilistic structure of the 
output current becomes very transparent. 

In Sec. |IV| we introduce the spectrum of the classical 
stochastic process which represents the output and we 
study the general properties of the spectra of such pro- 
cesses by proving characteristic bounds due to the Heis- 
enberg uncertainty principle. This bound is one of the 
evidences that the whole theory of continuous measure- 
ments is fully quantum, independently of the adopted 
formulation. 

As an application, in Sec. [V| we present the case of a 
two-level atom, which is measured in continuous time by 
detection of its fiuorescence light. The spectral analysis 
of the output can reveal the phenomenon of squeezing 
of the fluorescence light, a phenomenon related to the 
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uncertainty relations. We use this example also to il- 
lustrate the differences between homodyning and hetero- 
dyning and between the spectrum of the squeezing and 
the power spectrum. Finally we show how Mollow triplet 
appears in the power spectrum in the case of an intense 
stimulating laser. Section [Vl| contains our conclusions. 



we have 

bk{t)eif) = Mt)e{f). (3) 

It is a property of the Fock spaces the fact that the 
action on the coherent vectors uniquely determines a 
densely defined linear operator. 



II. THE QUANTUM STOCHASTIC 
SCHRODINGER EQUATION 

Quantum stochastic calculus is based on the use of 
some Bose fields playing the role of non-commuting 
noises and satisfying the canonical commutation relations 
(CCR) with a Dirac delta in time ([l]). This calculus and 
the unitary dynamics based on it were developed by Hud- 
son and Parthasarathy [TH], while Bose fields with delta- 
commutations in time were already found by Yuen and 
Shapiro ^29) in their study of the quasi-monochromatic 
paraxial approximation of the electromagnetic field. 



A. Quantum stochastic calculus and unitary 
dynamics 

Quantum stochastic calculus and the Hudson- 
Parthasarathy equation [T51 allow the evolution of a 
Markovian open quantum system, which we call system 
S", to be represented as a unitary evolution for system S 
interacting with some quantum fields. For a short review 
see [m Sec. 2] or Sees. 11.1, 11.2]; for a discussion 
of the physical approximations see [1^1 [H] . 



1. Bose fields 

Let us start by introducing the formal fields bk {t) , 6^ (t) 
satisfying the CCR 



b,{s),bl{t) ^S,k6{t-s), [6,(s),5fc(i)] =0. (1) 



In this paper we consider only the representation of the 
CCR ^ on the Fock space, the one characterized by the 
existence of the vacuum state. 

We denote by P = T {L"^ {R; C^)) the symmetric Fock 
space over the "one-particle space" L^(IR) — 
L2(|R; C^), and by e(/), / e ^^(IR; C''), the coherent vec- 
tors, whose components in the 0, 1, . . . , n, . . . particle 
spaces are 

e(/) := e-5 ll-^"' (l, /, (21)-^^' f ® /, . . . , {nl)-^^'f^\ . . .) 

(2) 

Note that e(0) represents the vacuum state and that 

(e(5)|e(/))=exp{-i 11/11^ -i||g|r + (.g|/>}. 

Let {zk, fc > 1} be the canonical basis in and for 
any / e L^iR^C^) let us set fkit) := {zk\f{t))c^. Then 



2. Factorization properties of the Fock space 

A symmetric Fock space P(3C) can be defined for every 
Hilbert space %; coherent vectors are defined always by 
([2]). When the one-particle space is given by a direct 
sum {% = %i © 3C2), one has the factorization property 

p(3Ci®3C2) = r(aCi)®p(3C2). 

In our set up, for every time interval A, let us de- 
note by T[A\ = P(L^(^;C'^)) the symmetric Fock space 
over L^(A;C''); in particular, we have P = P[IR]. Then, 
for any s < i, we have L'^{R]{:'^) = L'^{{-oo,s)]£,'^) ® 
L^{{s,t)-X'^) ® L^{{t,+oo)-X'^) and 

P[R] = p[(-oo,s)] ®P[(s,t)] «)P[(t,+oo)]. (4) 

Moreover, each space T[A\ can be identified with a sub- 
space of the full Fock space P[1R] by taking the tensor 
product of a generic vector in T[A\ with the vacuum of 
P[1R \ A]. Then, for every / e L'^{R]C'^), we have the 
identification 

e(/U)eP[A]K^e(U/)eP[R]. 

We are denoting by 1a(0 the indicator function of the 
set A and by /|a the restriction of the function / to the 
set A. With an abuse of notation we write 

e(/) = e(l(-oo,s)/) ® e(l(s,t)/) ® e(l(t^+oo)/) • 

In particular, e(l(3^t)/) can represent a vector in P[IR] 
or in P[(s,t)] and we have the identification e(l(s_t)/) = 
e(0)(8)e(l(,,t)/) ®e(0). 



3. Temporal modes and Weyl operators. 

The free evolution of the bose fields is represented by 
the left shift in P 

Qte{f)^e{6tf), {ej){s)^f{s + t). 

Then, the action of the shift on the fields is given by 

ej bkis)et^bkis + t). 

By using the energy representation of the fields, 

1 



bk{v) 



V 27r JiR 

we get <dl bkiy) 84 = e~"^*Sfc(i^), which shows that indeed 
Qt is the usual free dynamics. Therefore, the argument 
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t in the fields bk{t) has a double role: it is the time, 
because it appears in the evolution operator, and it is a 
field degree of freedom, the conjugate momentum of the 
free field energy. 

When the fields bk{t) represent the electromagnetic 
field, the index k can be used to denote the polariza- 
tion, the direction of propagation (discretized) or other 
spatial degrees of freedom 30J . 

If we take a function g E i^(IR) we can define the 
annihilation operator 

/ + 00 
g{t)bkit)dt. (5) 
-oo 

By Eq. (|3| , its action on the coherent vectors is given by 

/ + 00 
g{t)Ms)dseif) = {g\fk)mR)e{f). 
-oo 

If we take a complete orthonormal system g^, i = 
1,2,..., in i^(lR), we can define the annihilation oper- 
ators Ck{g^)- Together with their adjoint operators, they 
satisfy the usual CCR. We can say that the upper index 
i denotes the temporal modes, while the lower index k 
denotes the polarization/spatial modes. 

An important technical tool is represented by the Weyl 
operators VV(g), q € L^([R;C^), the unitary operators 
defined by: V/ S L'^{^-£'^). 

W(g)e(/) = cxp {iIm(/|(7)i2(|R.cd)} e(/ + q); 

this is nothing but the displacement operator for the field. 
By using the notation ([s]) we can write 

W(9)=cxp|^(cl(gfc)-h.c.)|, (6) 

while, by using the discrete modes introduced above, we 
have 

= exp 1^ ((.91%)l^(r)4(5*) - h.c.) | . 
By h.c. we denote the Hermitian conjugate operator. 



4- The Hudson-Parthasarathy equation 

Let M be the system space, the complex separable Hil- 
bert space associated to the quantum system S, and we 
take r, with its free evolution Qt, as the environment 
space. Now we want to construct the unitary evolution 
of the composite system on'K®V. 

By formally writing 

Bu{t) ^ f bk{s)ds, Blit)^ f\lis)ds, (7) 
Jo Jo 



we get the annihilation and creation processes, families 
of mutually adjoint operators, whose actions on the co- 
herent vectors are given by 

Bk{t)e{f)^ f h{s)dse{f), 
Jo 

{e{g)\Bl{t)eif)) = T ^ds (e(<?)|e(/)). 
Jo 

The overline denotes the complex conjugation. 

For t > 0, the annihilation and creation processes are 
adapted, in the sense that they factorizes, with respect 
to Q, as 

B';^\t) = l(_^,o) ® Bl^\t) (g) l(t,+oo), 
and they satisfy a variant of the CCR, namely 

[Bkit),Bj{s)]^5kitAs, (8) 
[Bkit),Biis)]^0, [Bl{t),Bj{s)]=0; 

t A s is the minimum between t and s and lyi is the 
identity operator on T[A]. 

By defining integrals of Ito type with respect to the in- 
crements of the quantum processes Bk, bJ,, it is possible 
to construct adapted operator processes on ?{ (g) F and 
to develop a quantum stochastic calculus, whose rules 
are summarized, at a heuristic level, by the quantum Ito 
table 

dBk (t) dBj (t) = Ski dt, dBl (t) dBi{t)^0, (9a) 
dBk {t) dBi{t)=0, dBl (t) dBj {t)=0, (9b) 
dBlit) dt = 0, dBk{t) dt = 0, (dt)2 ^ 0. (9c) 

Let Hq, Rk, k,l = 1,. . . ,d, be bounded operators on 
IK such that = Hq. We set also 

k 

Then, the quantum stochastic differential equation 
(quantum stochastic Schrodinger equation or Hudson- 
Parthasarathy equation) [THl HO] 

dUt = |^i?fcdB^.(t)-^i?^dBfc(i)+Xdt|c/t, (11) 

^ k k ^ 

with the initial condition C/q = 1, has a unique solution, 
which is a strongly continuous adapted family of unitary 
operators on J{(8)F, representing a system-field dynamics 
in the interaction picture with respect to the free field 
evolution [5T] . 

Then, for t > 0, the dynamics in the Schrodinger pic- 
ture is e~'^'^°'^* = QtUt, a strongly continuous unitary 
group whose Hamiltonian H^qt is a singular perturba- 
tion of the unbounded generator of Qt [321 133] ■ Roughly 
speaking, the system S is hit by a flow of bosons which 
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can have only a singular interaction with S; then, they 
are carried away by their free dynamics and never come 



back. The physical approximations involved in Eq. (11) 



are discussed in [51] Sec. 11.1.1]. Note that the interac- 
tion picture with respect to the free field dynamics co- 
incides with the Schrodinger picture when only reduced 
system states and observables are considered. 



B. The reduced dynamics of the system 

The states of a quantum system are represented by 
statistical operators, positive trace-class operators with 
trace one; let us denote by S(IH) the set of statistical 
operators on JC. For every composed state E in §(!H(g)r), 
the partial trace Trr (resp. Trj<;) with respect to the field 
(resp. system) Hilbert space gives the reduced system 
(resp. field) state Trr S in S(J£) (resp. Trjc E in S(r)). 

1. The initial state and the reduced states 

As initial state of the composed system "S* plus fields" 
we take gr{f) & S(iH (g) F), where p e §(?{) is gen- 
eric and gr{f) is a coherent state, gr{f) '■= |e(/))(e(/)|. 
Then, the system-field state at time t, in the field inter- 
action picture, is 



(12) 



We introduce also the reduced system state and the re- 
duced field state: 

77t:-Trr{E/(i)}, n/(i) Trjc {Ey(i)} • (13) 



2. The master equation 

One of the main properties of the Hudson- 
Parthasarathy equation is that, with the initial state in- 
troduced above, the reduced dynamics of system S ex- 
actly obeys a quantum master equation [THl HOI [12] . In- 
deed, we get 



(14) 



where the Liouville operator C{t) turns out to be given 
by 

C{t)[p]^^\[H^ + Hf{t),p] 

+ (RkpRl - \RlRkP - \pRiR,) , (15) 



thanks to the properties of the interaction and of the 
initial state of the environment. 

It is useful to introduce also the evolution operator 
from s to t by 



-T(t,s) = /:(t)oT(i,s), T(s,s) = l. (17) 
With this notation we have rjt = T(t, 0)[p]. 



III. CONTINUOUS MONITORING 

The connections among quantum stochastic calcu- 
lus, quantum Langevin equations and input and out- 
put fields were developed by Gardiner and Collet in 
[19) . Then, in [5 these notions were connected to 



the unitary evolution (111 and to continuous measure- 
ments. Indeed, another fundamental property of the 
Hudson-Parthasarathy equation is that it allows for a 
fully quantum description of a continuous measurement 
of the system S: the measurement is obtained by a de- 
tection of the bosons that have already interacted with 
S. Of course such a measurement acquires information 
on both S and the detected bosons. 



A. Input and output fields 

Let us call "input fields" the fields Bk{t), B\^{t),... 
when they are considered as operators in interaction pic- 
ture at time t, with respect to Ot, and let us call "output 
fields" the same fields in the Heisenberg picture: 



Bt\t) := UlBk{t)Ut 



(18) 



and a similar definition for B°™ (t) 



By the properties 
of the Fock space F and of the unitary operators Ut , it is 
possible to prove that 



BT\t) = UlBkitpT , 



VT> t. 



(19) 



This equation is of fundamental importance and it im- 
mediately implies that the output fields satisfy the same 
commutation rules of the input fields, for instance the 
CCR (|8|: the output fields remain Bose free fields. By 
applying the formal rules of QSC ^ , we can express the 
output fields as the quantum stochastic integrals [S| 



BT\t)^Bk{t)+ fulRkU^ds- 
Jo 



(20) 



given by the adjoint expression. 



H}{t) :=i^/fe(0i?fe-i^/fc(04- 



(16) 



B. The field observables 



Therefore, S is an open system, as it interacts with the 
fields in F, and its evolution turns out to be Markovian 



The key point of the theory of continuous measure- 
ments is to consider field observables represented by time 



5 



dependent, commuting selfadjoint operators in the Heis- 
enberg picture [H O [22] . Being commuting at different 
times, these observables represent outputs produced at 
different times which can be obtained in the same exper- 
iment. Here we present a very speciai case of observables, 
some field quadratures. Let us start by introducing the 
selfadjoint operators 



Q{t;^,h) 



h{s) dBlis) +h.c., t>0; (21) 



the phase i9 £ (— tt, tt] and the function h, with \h{t)\ = 
1, are fixed. The operators (21) have to be interpreted 
as linear combinations of the formal increments dB\{s), 
dBi(s) which represent field operators in the interaction 
picture. The corresponding operators in the Heisenberg 
picture are 

g°"*(t;i?,/i) : = UlQ{t-'d,h)Ut 

= UlQ{t;'d,h)UT, VT > t, (22) 



where the second equality follows from Eq. (19 1. These 
"output" quadratures are our observables. 

Each quadrature (5°"'(t; d, h) is observed at time t and 
it regards those bosons in "field 1" which have eventually 
interacted with S between time and time t, so it can 
be interpreted as an indirect measurement performed on 
the system S. 



By using CCRs, one can check that the operators (21 ) 
commute; \Q{t]'3,h),Q{s;'d, h)] = 0. The important 
point is that, thanks to Eq. ( |22| ), these operators com- 
mute for different times also in the Heisenberg picture. 
Therefore, the observables Q°"'(t; i?, /i), t > 0, can be 
jointly measured for every interaction (11). The output 
is a (random) number at every time i, that is a signal de- 
pending on time, a stochastic process, which is the result 
of a continuous indirect monitoring of the system S. Its 
probability distribution is given by the usual postulates 
of quantum mechanics trough the joint diagonalization 
of the operators Q°^^ {t; , h) . Actually, always thanks 
to Eq. (22 1, it will be enough to jointly diagonalize the 



operators Q{t;-d,h). 



Let us stress that quadratures of type (21) with dif- 



ferent phases and h functions represent incompatible ob- 
servables, because they do not commute but satisfy 

[Qit;^,h),Qis-ip,g)]^2i dr Im (e't''-'^)^^ ffW) • 

Note that for g = ft, we get 

[Q{t; 1?, h), Q{s; ^, h)] = 2i (t A s) sin (^ - y') , (23) 

and for tp — -d they commute as anticipated. 

When "field 1" represents the electromagnetic field, a 
physical realization of a measurement of the observables 
( [22| ) is implemented by what is called balanced hctero- 
dyne/homodyne detection [HI Sec. 8.4.4]. The 

light emitted by the system in the "channel 1" interferes 



with an intense laser beam represented by h, the local os- 
cillator. The mathematical description of the apparatus 
is given in |2U Sec. 3.5]. 

Let us note that the operator Q°"*(t; d, h) involves the 
whole time interval [0, t] and has to be interpreted as 
cumulated output. The instantaneous output current 
is represented by its formal time derivative /°"*(t) :— 
Q°'''it;i),h). From (0, ([l9|, ^ we get 



UjRiUt ] +h.c. 



(24) 



C. The stochastic representation 



The commuting selfadjoint operators (21) have a joint 
projection valued measure (pvm) E^, which gives the 
probability distribution for the output of the continu- 
ous measurement. Moreover, via the partial trace on 
the fields, gives also the instruments describing the 
transformations of S from time to an arbitrary time 
t, conditioned on the information acquired up to time 
t. Furthermore, via joint diagonalization and condition- 
ing, the pvm E^ even gives the stochastic evolution of 
the conditional state pt (or a posteriori state), the state 
of S at time t given the observed signal from time to 
time t. This evolution turns out to satisfy a stochastic 
differential equation (SSE or SME), with classical driv- 
ing noises. The introduction of such stochastic evolution 
equations for the conditional state was an achievement 
of the quantum filtering theory [71 137H40| . 

The passage from the formulation with quantum fields 
and Hudson-Parthasarathy equation to the one based on 
classical stochastic differential equations can be done by 
different techniques. The technique based on the use of 
isomorphisms between the Fock space and the Wiener 
space is very powerful and clear; here we present a variant 
of the construction given in [12 . 

Let us note that the observation we consider is not 
complete, because it regards only field 1 and involves only 
positive times. To make unique the isomorphism which 
diagonalizes the self-adjoint operators (21 1, we need to 



add fictitious observations, involving quadratures of the 
fields 2,...,d too. So, we take a function £ e L°°{R; C*) 
such that 

14(01 = 1, Vt e R, Vfc, and £i{t) = e-'^h{t), Vt > 0; 
then, we introduce the field quadratures: for k = 1, . . . , d. 



Qk{t) 



£kis)dBl{s) +h.c. 



(25) 



We use this definition for positive and negative times by 
taking the convention — — J^" for a negative t. These 
quadratures form a complete set of compatible observ- 
ables on the Fock space F. Note that Qi{t) = Q{t;d,h). 
In the following subsection we jointly diagonalize all the 
observables ( 25 ) by introducing an explicit isomorphism 



between Fock and Wiener spaces. 
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Spectral representation on the Wiener space 



Fixed the functions £i,. . . ,£d, that is the field quad- 



The distribution of the output. Let us restrict 



ratures (25), we look for a probability space (fi, J, Q), 
a unitary operator J : T[L'^{R;C'^)] L2(rj, (the 
Hilbert space of the complex square integrable random 
variables on the given probability space), and a family of 
random variables Wk{t) on il such that 

{jQkm){ij) = Wkit;ij) (J*)(c^), (26) 

for all t, k, and for all in the domain of the selfadjoint 
operator Qk{t). This means that each Qk{t) is repres- 
ented in L^(ri, 3^,(5) as the multiplication operator by 
Wk{t). We can get such a joint diagonalization on the 
space of the canonical representation of the Wiener pro- 
cess; the canonical Wiener process is presented, for in- 
stance, in [m Sees. A.2.4, A.2.6]. 

Let r2 = Co(lR; R'') be the space of the continuous 
functions w : IR — )• R'' such that a;(0) = 0. We define 
the d-dimensional process W{t) : Vl — IR'*, t G R, by 
W{t,uj) = ijj{t) and we denote by 'J the smallest cr- 
algebra of subsets of Vl for which these functions W{t) 
are measurable: J — a{W{t) : i € R). Then, there 
exists a unique probability measure Q on the measurable 
space (rj, 7), the Wiener measure, such that the processes 
Wk{t), Wk{—t), t > 0, k = 1, . . . ,d are 2d independent 
standard Wiener processes. Moreover, for positive times 
we introduce the natural filtration (3^t)t>o of the process 
W: Jt = cr{W{s) : s € [0,t]). Let us also recall that, if 
4>,tp G L^(r2,S', Q), then their inner product is given by 
the Q-expectation Eq: 



(V'|0) = Eq[7A0]= / V(^)'/'(^)Q(dc^ 



Let J : r[L2(R; C*)] ^ ^^(^2, J, Q) be the linear oper- 
ator defined by: V.g G ^^(R; C''), 




X exp 

^"•^-<^ 

In particular wc have 

Je(0) = l, Je(l(o,t)/) eL2(J7,?,,Q). 

The operator J turns out to be an isomorphism and it 
reahzes the representation (26 1: J Qk{t) J^^ = Wk{t)^ 
i.e. the field quadratures are mapped into the operators 
"multiplication by the Wiener processes". Because the 
isomorphism J jointly diagonalizes all the observables 
(251, then their joint pvm on the Fock space is J~^1^J, 



now to the observed quadrature (21 1; the tr-algebra 
Soo = <^{Wi{t) : Q < t < -foo) is the space of all the 
events regarding our observables Q{t; h), t > 0. Then, 
the joint pvm of the observed quadratures is defined 
on the measurable space (f2, Soo) by 



(28) 



Finally, we get the distribution of the output. By setting 
9t = cr{Wi{s) : s e [0,f]), then, St C Soo, is the space of 
the observed events up to time t, associated to the ob- 
servables Q{s; h) for times from to t, and, according 
to the usual rules of quantum mechanics, the probabilit- 
ies of such events are given by 

P^t(G) = Tr{(lIjt®£;^(G))S/(i)}, VG G St, Vi > 0. 

(29) 

Note that, when the field state is the vacuum and there 
is no interaction between system S and the fields, this 
probability reduces to Pp;t'(G) = {e{Q)\E^{G) e{Q)) = 
Eq[lG] = Q(G). This means that in this case the quad- 
ratures (pTl) are distributed as a standard Wiener process. 



Let us stress that the pvm (281 depends on the para- 



meters •& and h defining the quadrature (21 ); these para- 



meters are contained in the definition of the isomorphism 
J (27). On the contrary, the choices of the trajectory 



space (the measurable space (51, Soo)) and of the pro- 
cess Wi are independent of d and h. With respect to 
the time dependence, the physical probabilities (29 1 are 
consistent i.e. 



< s < t, G e S. 



=p:f (G) = p;^:,''(G). 



(30) 



This result is due to the factorization property (|4| of the 
Fock space and to the localization properties of Ut 1211 
Theor. 2.3], which imply u} {\.k® E'^{G))Ut = Ul{\^® 
E'^{G))Us for < s < i and G G Ss, cf. Eq. ([22]). 

A more detailed study of the statistical properties of 
the output needs the introduction of the characteristic 



operator (Sec. HID) 



2. The instruments 

The observation of the emitted field can be interpreted 
as an indirect measurement on the system S and this 
is formalized by the concept of instrument [2 H] . The 
family of instruments Tt, t > 0, describing our measure 
is defined by: VG £ St, Vr e S(J{), 

It(G)[T] = Trr {(Ijt ® E^{G))Ut{T <E> gr{f))u}} . 

(31) 

For T — p, the initial system state, Eq. (31) gives the 



non-normalized state of S at time t conditioned on the 
information that the values of the signal in the time in- 
terval from to t were in G. Of course we have 



Tr^{It{G)[p]} 



(G), 



(32) 
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while the normaHzed conditioned state is given by latter instrument has the simple a posteriori state (38 1, 



Tt{G)[p] divided by its trace (32) 
Let us remark that 



(33) 



so that the system reduced state at time t in the case 
of no observation (ry^) coincides with the so called a pri- 
ori state (It(r2)[p]), that is the system state at time t 
in the case of observation performed but not taken into 
account. This is in agreement with our rough picture of 
the measurement process: we observe fields which have 
already interacted with system S and which will never 
interact again with it. This means that we acquire in- 
formation on S", as we have It{G)[p] ^ P^'f (G)77(, but we 
do not add any perturbation on its evolution as we have 

a. The a posteriori states. Now we want to intro- 
duce pt^ the state of S at time t conditioned on the 
whole information supplied by our indirect measurement 
between time and time t, that is by the signal produced 
by the measurement of Q{s] d, h) for s E [0, t]. Therefore, 
pt has to be a random state depending on the output 
Wi{s), < s < t, that is a random state measurable 
with respect to St; in other terms, we have the func- 
tional dependence pt(w) = pt(a;i(s), < s < t). Such 
a state is called a posteriori state and it is determined 
by the initial state p and by the instrument It: it is the 
unique Sf-measurable random state such that 



MG)[p] = 



Pt(c.)p!.f(dc.), 



VGe St 



(34) 



The definition of a posteriori state is not linked only to 
measurements in continuous time, but it has been intro- 
duced for a generic instrument |41| . 

As we have a reference probability Q on the out- 
put space (f2,St) we can equivalently look for the 
non-normalized a posteriori state at, the unique St- 
measurable random positive operator such that 



It{G)[p] = / <Jt{uj)€i{dij), VG€ St 



(35) 



Then, Trjat} is the probability density of P^'f with re- 
spect to Q and we have pt — df/ Trjcrt}. 

The non-normalized a posteriori state at can be com- 
puted by using the spectral representation ( 26 1 of the op- 



erators Qk and its evolution can be obtained by passing 
through the SSE. 



3. The stochastic Schrodinger equation 

In order to compute the a posteriori state of our in- 
strument Xt (31), it is convenient to pass through two 
fictitious instruments: J7f, associated to a complete set 
of compatible observables in T, and J7f, associated to a 
complete set of compatible observables in r[(0,i)]. This 



whose evolution is given by the SSE (40 1 



First of all, let us imagine, in the Heisenberg picture, 
that in the time interval [0, t] we measure all the quad- 
ratures Qr'(s) = UjQkis)Us, k = l,...,d, s £ [0,t], 
and moreover we conclude the measure by observing at 
time t also the field observables Qk{u;t) — Ut Qk{u)Ut, 
k — 1, . . . , d and u < or u > t. This is a family of 



commuting observables, thanks to (23) and to (19), that 
U}Qk{s)Ut. Then, the instrument jTt 



implies Q^"*(s) — ^t '^k\^)^t- inen, me iiisu uiiieiii Jt 

associated to this fictitious measurement is given by an 
expression analogous to (31). If the system initial state 
is pure, p = |r)(r|, r G 3^77^1 



1, then VF € J, 



Jt{F)[\r){r\]=Tvr {{t-K® J-^lFJ)\^t){^t\} , (36) 
where 

*t = t/t(r®e(/)). 



To get the expression ( 36 1 we have changed picture and 



cycled Ut against the state; then we have used the iso- 
morphism J (27) that diagonalizes all the quadratures 
Qk{u) 

The isomorphism J^^ does not involve the space J{ 
and it can be cycled after (5't|; in this way we get 

Jt{F)[\r){r\]^ f |¥>t(^))(^t(a.)|Q(da.), (37) 
Jf 

where ipt is the random IK-vector 

ipt = J^t = JUt{r®e{f)y 



By comparing Eq. (37) with Eq. (35), we get that the 



non-normalized a posteriori state af associated to the 
instrument and to the premeasurement system state 
p = \r){r\ is 



erf 



(uj) = \(ptiuj)){(pt{(^)\- 



By construction af is a random positive trace-class op- 
erator, which is ^F-measurable. 

Suppose now that we measure only the quadratures 
Qr'(s)> < s < k = l,...,d; note that this set 
of compatible observables is complete in r[(0,t)], not in 
r[R]. With respect to the previous case, we simply have 
to drop some commuting observables and thus the new 
instrument J7t is just the restriction of J^t to the u-algebra 
ff't C 9^ and, therefore, the new a posteriori state af is 
the conditional expectation 



which is an ff't -measurable random positive operator. 
Thanks to the properties of the Hudson-Parthasarathy 
equation and to the choice of the observed quadratures 
(local in (0, t) and no k neglected) the a posteriori state 
is still almost surely pure: 



afioj) = \M^)){M^)\, ^t = J(7t(r®e(/l(o,t))). 



(38) 
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a. The linear SSE. It is the (^Fi)-adapted stochastic 
process (f)t that satisfies the hnear SSE and we show now 
how to get it. 

By introducing the Weyl operators Wt ■= W(/l(o,t)) 
we can write e(/l(o.f)) = yVfe(O) and 



(f>t = JUtWt{r(g)e{0)). 



(39) 



By using the definition of the Weyl operators and the 
quantum stochastic calculus it is easy to check that Wt 
satisfies a quantum SSE ( 11 ) with = fk{t) and Hq = 



0. Then, the quantum stochastic differential of UtWt is 
given by 

d {UtWt) = {{dUt) Wt + Ut dWt + (dUt) dWt) ■ 

This can be computed by using the quantum Ito table ^ 
and exploiting that operators localized in disjoint time 
intervals commute and that the differentials dB'j^\t) are 
localized in (t,t + dt) with respect to the factorization Q 
of the Fock space; in particular this means that dBk[t) 
and Ut commute and the same holds for dBk{t) and Wf 
The result is 



d (UtWt) = Udt + J2 [{Rk + fkit))dBiit) - (4 + Mt)) dB.it) - (i \Mt)f + hmi) dt 



UtWt 



Now, we use this result to compute the differential of 0f ( 39 ) . The key point is that 



{dBk{t))Ut Wt (r (g> e(0)) = Ut Wt dB^lt) (r ® e(0)) = 0, 
so that we can change the coefficient of dBk{t) as we wish. By this and the fact that = 1, we can write 



= J I Kdt 



51 [(^fe + f kit)) hit) dQkit) -[l\fkit)\^ + fkimi) dt\ \ut Wt{r ^ e(0)). 



Finally, by Eqs. (10), (26), and (39) we obtain the linear SSE: 



d<j>t = ik{t){Rk + fk{t))dWk{t) - ^ (4 + fkit)) {Rk + fk{t))dt -i Ho+'-J2 {h{t)Rk - fk{t)Rl 



dt 



(40) 



Therefore, the a posteriori evolution of system 5* under 
the continuous measurement Xt is a stochastic evolution 
mapping pure states into pure states. In particular, the 
evolution of the non- normalized pure state 4>t is given by 
the linear SSE (40) and it is Markovian and depends on 
the interaction ( 11 ) between S and F, on the field initial 



4-. The stochastic master equation 



state grif) and on the observed quadratures Q1 (s). 



It is now possible to show that ||0x(a;)||r^ Q(dw) defines 
a new probability on {Q^S't) and that, under this new 
probabihty, the process ■— 0t/ll0t||jti t ^ i^^T], sat- 
isfies a nonlinear stochastic differential equation. This 
last equation is the nonlinear SSE, which is the start- 
ing point for useful numerical simulation. A key point 
in the change of probability is the fact that ||0t||^ is a 
Q-martingale and that the so called Girsanov transform- 
ation can be invoked. For the theory of the linear and 
nonlinear SSE we refer to jTHl Sec. 2]. 



By Ito calculus, from the SSE (40 ) we get the stochastic 
equation satisfied by the random operator ([38]): 



daf ^Cit)[af] dt + Y,{h{t){Rk + fk{t))af 

k 

+ <jf 4(t) (rI + W)) ]dWk{t), (41) 



where C{t) is the Liouville operator (15). 

We can now get rid of the hypothesis of a pure ini- 
tial state and prove that Eq. (41) gives the a posteriori 



evolution for a generic system initial state 



1, > 0, ^p£ = l. 



I 

Indeed, if we set 

(l>l^JUtWt{r,®e{Q)), 
then, by linearity the process 

(^f {^) ^^viWt{^)){^l{^)\ 
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is adapted, satisfies Eq. (41), and gives the a posteriori 
state of Jt, 



Jt[F)[p] 



a/HQ(dL^), 



(42) 



Finally, we consider our instrument Zj, which is the 
restriction of Jt to St, so that Eq. ( [35| , defining the non- 
normalized a posteriori states, holds with 



(Jt = E( 



(43) 



Let us recal that ii{t) ~ e h{t)^ that the (Q-mean of any 
Wk[s) is zero, that St is generated by Wi and that the 
other components of the Wiener process are independent 
from the first one. Then, by applying the conditional 



expectation with respect to St to (41), we obtain the 



linear SME for the non-normalized a posteriori states of 



C{t)[at]At+[c'''h{t){R^+h{t)) 



+ crt e" 



-'^h{t)(R\+h{t))]AW^{t). (44) 



As already said at the end of Sec. |III C 2[ the quant- 
ity Trjt {(Tt(w)} is the density of the physical probability 
with respect to Q. Indeed, from Eqs. (35) and (32) we 
get 



p,t 



(G) 



TrM{tTt(w)}Q(dw), VGeSt. (45) 



G 



By taking the trace of the SME (44), it is possible to 
show that Trj^- {at} is a Q-martingale. In the stochastic 
formulation it is just this fact that implies the consistency 
property ( 30 ) , which we already encountered in the Fock 



space formulation. 

As already seen, if we define pt = (Jt/ Trj^jcTf }, we get 
the a posteriori state for the instrument Tt and the pre- 
measurement state p. By ( 33 ) , the a posteriori states are 



related to the system reduced state by 



pt{uj) P^,t (dw) = / atiio) Q(dc.) = ,7t. (46) 



It is also possible to prove that, under the physical 
probability, pt satisfies a nonlinear SME. Moreover, by 
studying the stochastic differential of Trjf {at} and by 
using Girsanov theorem, it is possible to prove the fol- 
lowing result. 

a. The noise. Under the physical probability Pp'^ 
the process 



Wi{t):^Wi{t)-2^e / e'''h{s)TY^ {{Ri + h{s)) p,] As 

JO 

is a standard Wiener process for t £ [0, T]. 

In other terms we can say that, under the physical 
probability, the instantaneous output I{t) = Wi{t) is 



the sum of a white noise AWi (t) / At plus a regular signal 
2ReTr5f {{Ri + fi{t))pt}- White noise and signal turn 
out to be correlated in general. Let us stress that this 
result on the structure of the output is a byproduct of 
the stochastic representation of the continuous measure- 
ments. From this representation and Eq. (46) the mean 



value of the observed quadrature at time t is 



Tr{(ljt®Q(t;z9,/i))S;(0} = / Wi{t; u;) P;f (Au;) 

Jn 

= 2Re f e"'Ms)TrM{(i?i+/i(s))77s}ds. (47) 

JQ 

The full theory of the linear and nonlinear SME's and 
their relations with the physical probability and the ref- 
erence probability Q are presented in [TBI Sees. 3 and 5]. 



D. Characteristic functional and moments 

In the study of stochastic processes it is often useful 
to have explicit formulae for the moments, for instance 
for the second order moments, which determine the spec- 
trum of the process; see Sec. |IV A] In the case of our out- 
put, the mean function is given by Eq. (47); to get the 



higher moments it is useful to introduce the character- 
istic functional, which is the functional Fourier transform 
of the probability distribution of the process. 

Let us denote by E^'f the expectation with respect 

to the physical probability Pp'f. By recalling that the 
output is represented by Wi , the characteristic functional 
up to time t > is 



$t(fc;z9,;i) = E 



P.t 



exp 



1 / k{s)AWi{s) 





(48) 



the argument k is any real test function in L°°(R+). 

By functional differentiation with respect to the test 
function one gets all the moments of the process, as 
done in Sec. |III D 2| Moreover, when the characteristic 
functional is given, one gets the probabilities by anti- 
Fourier transform. For instance, the finite-dimensional 
probability densities of the increments Wi{ti) — Wi{to), 
Wiit2) - Wiiti), VFi(t„) - Wi{tn-i), with 0<to< 
ti < ■ ■ ■ < tn < t, are given by 

where we have introduced the test function k{s) — 

E"=i i(t,-i,t,)(s)«j- 



1. Characteristic operators 

As our measurement is an indirect observation of S 
performed by a direct observation of F, the character- 



istic functional ( 48 ) can be expressed either in terms 
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of the system S only or in terms of the fields only. 
First, we want to express it in terms of the quantum 



observables (21). We introduce the characteristic oper- 
ator ^t{k; t?, h), the Fourier transform of the pvm [51] 
Sec. 11.4.2], 1211 Sec. 3.2]: 



(49) 

By using the representation ( 28 ) of the pvm, the cor- 
respondence between Q and Wi (26), and the definition 
of Q (pTl), we get 



2. The output moments 

By functional differentiation of the characteristic func- 
tional, we get all the moments of the classical output 
process. Let us introduce the formal time derivatives 
I{t) = Wi{t) and i{t) = Q{t;§,h); from ^ and ^ 
we obtain immediately the expressions of mean function 
and autocorrelation function: 

E^;^[/(t)] = Tr{Q(t;z?,/i)S;(T)} = Trr {/(<)%(T)} 

= 2Re(e'''/i(t)Trr{6i(i)n/(r)}) , (53a) 



$t(fc;?9,/i) = J-iexp|iy k{s)dWi{s)}J 
= cxp|i J k{s) dQ{s;'&,h) 
= cxp|ic-'''^ k{s)h{s)dBl{s) -h.c. 



(50) 



By comparing this expression with Eq. (|6|, we see that 
the characteristic operator is the unitary Weyl operator 

$t{k;^,h) = Wiq), qj{s) = S,iie'''>k{s)h{s)l[o,t]{s). 



Then, by using the expression (29) of the physical 



probability, the characteristic functional can be written 
as 



:(fc;i?,/i) =Tr{$t(fc;i?,/i)E/(t)} 



Tr < exp 



i / /c(s)dg™'(s;z9,/i) 



£'r(/) 



Trr{$t(fc;t?,/i)n/(i)} 



(51) 



The last step is due to the fact that $t(fc; -d, h) depends 
only on field operators; recall that S/(t) is the state of 
the total system S plus fields, while 11/ (t) is the reduced 
state of the fields ( 13 ). 



Finally, let us define the reduced characteristic operator 
Qt as the functional Fourier transform of the instrument 



^ [aniiniisa: 

gt{k;d,h) = / exp<; i 



k{s) dWiis)j{Lo)'iIt{duj). 

(52) 

It can be shown that Qt satisfies a closed differential equa- 
tion, a kind of modification of the master equation 0]. 
Then, by the representation ( 32 ) of the physical probab- 



ilities, we get a further expression of the characteristic 
functional: 

<S>t{k;d,h)^TT:K{gtik;^,h)[p]}. 



5{t - s) + 2Re(his)Trr[{h{t)b\it) 
e^'^W)bi{t))h{s)Uf{T)}), 



(53b) 



where T > t, T > s. Analogous formulae hold for 
higher moments. Let us note that the order of the op- 
erators I{t) and /(s) in (53b) does not matter, because 
they commute. Moreover, the moments of the classical 
process I{t) are expressed in terms of quantum means 
and quantum correlations of the fields [211 P- 165 and 

Sec. 11.3.2]: Trr {6i(i) n/(r)}, Trr {6l(i)&i(s) n/(r)}, 

Trr {6i(i)6i(s) n/(r)}, and the complex conjugated ex- 
pressions. The fields are all in normal order because we 
put in evidence the delta term coming out from a com- 
mutator. 

By studying the properties of the reduced character- 
istic operator (52 ), it is possible to prove that all the mo- 



ments of our classical output can be expressed by means 
of quantities concerning only system S [22j . For the mean 
and autocorrelation functions the final result is [55J Sees. 
3.3, 3.5] 



E;;^[/(t)]=2Re(Trjt{Z(0r;J), 



E;;^[/(t)/(s)] - Sit - s) + 2Re(Trjt{Z(t2) 



(54a) 



X T(t2,ti) [Ziti)rjt,+Vt^Ziti)^]}), 

(54b) 

where t2 — t\/ s, ti = t A s and 

z{t) -.^c'^^W) (Ri + fiit))- 

These expressions are more useful for computations, 
while the expressions ([53]) are better suited for theor- 
etical considerations; cf. [21[ Sec. 5.4.6]. 

IV. THE SPECTRUM OF THE OUTPUT 

Inside the theory of continuous measurements, the out- 
put of the measurement is a classical stochastic process. 



11 



whose distribution is determined by quantum mechan- 
ics; so, the spectrum of the output can be introduced by 
using the classical definition of spectrum of a stochastic 
process 



A. The spectrum of a stationary process 

In the classical theory of stochastic processes, the spec- 
trum is related to the Fourier transform of the autocorrel- 
ation function [l^ . Let y be a stationary real stochastic 
process with finite moments; then, the mean is inde- 
pendent of time E[Y{t)] = E[Y{0)] =: my^^t & R, and 
the second moment is invariant under time translations: 
Vi, s e R, 

E[r(t)y(s)] = E[Y{t - s)r(O)] =: Rvit - s). (55) 

The function Ryit), t € R, is called the autocor- 
relation function of the process. Obviously, we have 
Gov [Y{t), Y{s)] = Ryit - s) - m^. 

The spectrum of the stationary stochastic process Y is 
the Fourier transform of its autocorrelation function: 



Sy{fl) : = 



+CX3 



e'^''Ry{t) dt. 



(56) 



This formula has to be intended in the sense of distribu- 
tions. If Gov [Y{t),Y{0)] e L^{R), we can write 



Sy{ii) := 2'KmyS{ii) + 



+ 00 



a'*'* Gov[y(t),r(0)]dt. 



(57) 

By the properties of the covariance, the function 
Cov [Y{t),Y{0)] is positive definite and, by the prop- 
erties of positive definite functions, this implies 
e'^* Gov [Y{t), Y{0)] dt > 0; then, also Syifi) > 0. 
By using the stationarity and some tricks on multiple 
integrals, one can check that an alternative expression of 
the spectrum is 



T— !.+oo 1 



e'''*r(i) dt 



(58) 



The advantage now is that the positivity of the spectrum 
appears explicitly and only positive times are involved 



which are stationary only in some asymptotic sense and 
to singular processes as our I{t). 



B. The spectrum of the output in a finite time 
horizon 

Let us consider our output I{t) = dWi{t)/dt under the 
physical probability Pp'-^- We call "spectrum up to time 
T" of I{t) the quantity 



3"^* dl^i(t) 



(59) 



Note that the spectrum is an even function of /z: 
S't(m;i?) = 5'r(-M;^^)- When the limit T +oo exists, 
we can speak of spectrum of the output, but this existence 
depends on the specific properties of the concrete model. 

By writing the second moment defining the spectrum 
as the square of the mean plus the variance, the spectrum 
splits in an elastic or coherent part and in an inelastic or 
incoherent one: 



2 



^5J(/x;z9): 



1 

f 



1 i).h 



— Var'''™ 



3'^* dWi{t) 



cos ^t dWi (t) 



sin fit dWi{t) 



(60a) 
(60b) 



(60c) 



Let us note that 

s-^ifi; = 55J(-A^; ^), sr\i^\ ^) = 

(61) 

In Sec. |IIID| we have seen two ways of expressing the 
output moments, by means of field operators or by means 
of system operators. 



By using the expression (53b) for the autocorrelation 



function of the output, we obtain 

S't(m; + ^ f dt [ ds e'^('-") Re(h(s) 

X Trr { (h{t)b\ it) + e^'"" h{t) bi {t)) h {s) B; (T) } 



(62) 



Let us note that Eq. (62) expresses the spectrum as a 



Fourier transform (in a finite time interval) of a normal 
ordered quantum correlation function of the field; cf. jfSl 
Sec. 9.3.2]. 



By using the expressions ( 54 ) for the first two moments 



we get the spectrum in a form which involves only system 
operators: 

rr. 2 



Expression ( 58 ) can be generalized also to processes (fi; 



4 
T 



1 

f 



3'''*Trjt{(Z(<) + Z(0t)r,,}dt 



Re 



h{t){TY^{R,7Jt} + f,{t)) 



dt , 
(63a) 



S'^"°'(Ai;i9) = 1 + I; / dt / ds cos fi{t - s) 
^ Jo Jo 



X Tr 



:k { {Z{t) + Z(i)t) T(t, s) [Z{s)rfs + VsZ{s)^] } , 

(63b) 



Z{t) = Z(t)-TrM {Z{t)i^t} = e'" /i(i)(i?i -Trj^ 



12 



C. Properties of the spectrum and the Heisenberg 
uncertainty relations 



Equations ( 63 ) give the spectrum in terms of the re- 



duced description of system S (the fields are traced out); 
this is useful for concrete computations. But the general 
properties of the spectrum are more easily obtained by 
working with the fields; so, here we trace out system S 
and we start from expression ( 62 ) . 



Let us define the field operators 



1 



e'''* dQ{t;^,h), 



(64a) 



Qrifi; ^) = QT(/i; i?) - Trr {Uf{T)QT{fi; i9)} ; (64b) 

the local oscillator wave h is fixed. Let us stress that 
Qt{p-', '&) commutes with its adjoint and that Qxin', i?)^ — 
Qt{—IJ'\'&)- By using Eqs. (51) and (53) and taking first 
the trace over iK, we get 

ST{^i;'^) = TrT{Ilf{T)QT{^i■,'^yQT{^i■,■^)]>Q, (65a) 

S-^(p-d) = |Trr{ny(r)QT(A*;i?)}|' > 0, (65b) 

::rinel / 



S'^''\^i■ ^) = Trr {%(T)QT(/i; i9)^Qt(a*; i?)} > 0. 



1. Spectrum and field modes 



To elaborate the previous expressions it is useful to 
introduce annihilation and creation operators for bosonic 
temporal modes, as in Sec. |II A3[ 

aT(/i) 4f /'^e'^*Mt)dSi(i) ^ c^{gl^). (66) 
vJ Jo 



Ht)l[o,T]it). 



(67) 



The operators aT(/i), ^'^(/i) are true bosonic modes, as 
they satisfy the CCR 

[aTi^i), arifi')] = [a^(Ai), a^(Ai')] = 0, (68a) 
[aTi^l),a^^i^l)] = ^. (68b) 

However, for finite T these modes are only approximately 
orthogonal, as we get 



[aTi^J.),a^p{^l')] 



i(M-A^')r 



for fi' ^ ^. (69) 



Then, from Eqs. ( |2l| ), (64a), ( |65a ) we have easily 



ifi; = e'\T{^l) + e-'\U-f^)^ (70) 



STi^l■, ^) = Trr{ (e-"'a^(-Ai) + e'%T(Ai)) 

X n/(r) (e-"'a^(M) + e"'aT(-/^))}. (71) 



Let us stress that only two field modes contribute to the 
spectrum for fj, ^ 0, and only one mode in the case of 
^ = 0. By using the CCR ( 68 ) we get the normal ordered 
version of ( 71 ): 



5t(m;i?) = i 

+ TTT^Tlf{T)(^a\.{fi)aT{fJ.) + a^(-^)aT(-M) 



+ e 



-2it5„t 



al{fi)al{-fi) + e2'''aT(-/i)aT(Ai)) }. (72) 



Note that Eq. ( 69 ) played no role in the normal ordering 
operation. 



By Eqs. (65b) and (70) we get for the elastic part of 



the spectrum 

S^\^l■,^})^\TT^{UfiT)aTi^^)} 

+ e-2i''Trr{%(r)4(-Ai)} 



(73) 



To obtain a similar expression also for the inelastic 
part, it is convenient to introduce the operators 

aT(/i) aT(/i) - Trr {nf{T)aT{fi)} , (74) 



which satisfy the same commutation relations mSh , ( 69 ) 



(65c) as the operators arifJ-) and their adjoint. Then, we get 
Qt(m; ^) = e'^&rifi) + e-"'4(-/x), (75) 



+ Trr|n/(r)(^4(A^)^T(M) + a\:{-^i)aT{~^) 

+ e-2"'4(^)4(-^) + e^'''~aT{-^lYaT{^i))]. (76) 



2. Spectra of complementary quadratures 
Let us consider two choices of the phase d: d and § ± 



7r/2. From Eq. (j^Sj) we get 

[Q[t- h),Q{s; ± tt/2, h)] = T2i {t A s) , 

which means that we are considering two incompatible 
field quadratures, measured by two different setups. For 
these quadratures we have the important bounds and re- 
lations given in the following theorem. 

Theorem 1. For every "i? and fj, we have the following 
relations: 
1 



^ 5?}(M;^) + 5?J(A^;^±f)) 
= |Trr {nf{T)aT{ti)}\^ + |Trr {Uf{T)aT{-fi)}\' 



(77a) 



-(sr'(M;^) + ^r'(M;^9±f)) = 1 



Tr 



r { n/ (T) (^al, {n)hT (a*) -t- 4 (-A*)) } , 



(77b) 
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+ Trr |n/(r) (^5.'l{^i)dT{p) - d\,{-fj.)aTi-ti)j | 

(78) 

Then, independently of the system state p, of the field 
state Qrif), of the function h and of the Hudson- 
Parthasarathy evolution U , the following two important 
hounds hold: 



Equations ( 77 1 express the independence from d of the 



sr\^Ji■,^)ST\^l■^±l)>l, 



(^r'(M;^9) + 5T (M;^±f)) >i- 



(79) 



(80) 



Proof. First of all from Eqs. (73), (76) we get Eqs. (77 1 



Then, the bound (80 1 comes immediately from Eq. (77b) 
The bound ( 79 1 is a trivial consequence of Eq. ( 78 



To prove the bound ( 78 ) , we write 



S'T'ilJi: ^) = Trr{ (e-"'4(-Ai) + e"'aT(/^) 



n/(r) (e-'''4(Ai) + e'''5T(-A^))}. 



The usual tricks to derive the Heisenberg-Scrodinger- 
Robertson uncertainty relations can be generalized also 
to non-selfadjoint operators [S]. For any choice of the 
state Q and of the operators Xi, X2 (with finite second 
moments with respect to g) the 2x2 matrix with elements 

Tr|Xi(3Xj| is positive definite and, in particular, its 

determinant is not negative. Then, we have 



TT{XigXl}TT{x2gXl} > \Tr {x^gX^} 

> Im Tv{x,gXl}\^ = ^\Tv{g (^X^X^ - Xlx^) } 
By taking e = n/(r), 

=e-'^4(-Ai)+e'''aT(Ai), 



X2 = Tie-'''4(-M) ± ie'''&T(M), 

we get 

> 1 + Trr{n/(T) (4(-Ai)aT(-M) - 4(m)«t(m)) } 
But we can change fi in —fi and we have also 



> 



+ Trr {n/(T) (a^ {p.) drip-) ~ 4(-m)«t(-m)) } 



The two inequalities together give the final result (78 1. 

□ 



arithmetic mean, in both cases of elastic and inelastic 
spectra. Equation (78) is a bound of Robertson type; 



such a bound does not depend on but it is still de- 
pendent on the initial state and on the dynamics. The 
Heisenberg-type relation (79) and the bound (80) are 



fully independent of the initial state and of the dynamics. 

One speaks of squeezed field or of the spectrum of 
squeezing [151 Sec. 9.3.2] if, at least in a region of the 
/i line, for some one has S™'^^{fi;d) < 1. If this 



happens, the bounds (79) and (80) say that necessar- 
ily S^r^°^{fi;'d + f ) > 1 in such a way that the product 
and the arithmetic mean are bigger than one. Note that, 
with our choice of the environment initial state, any pos- 
sible squeezing can be imputed to the interaction with 
S. Indeed, in the case of no interaction, the output 
Wi is a Wiener process plus a deterministic drift and 
(A^,^)^l. 



cincl 



V. HOMODYNING VERSUS HETERODYNING 
OF THE FLUORESCENCE LIGHT OF A 
TWO-LEVEL ATOM 

Let us take as system S a two-level atom, which means 
5{ = C^, Hq = ^ cr^; > is the resonance frequency 
of the atom. We denote by (t_ and cr-|_ the lowering and 
rising operators and by ax — <y- + cr+, cry = i(tT_ — tT+), 
az = cr^a^ — (T_(T_|_ the Pauli matrices; we set also 



+ c-'%+, 



The atom can absorb and emit light and it is stimulated 
by a laser; some thermal environment can be present too. 
The quantum fields F model the whole environment. 

a. The absorption/emission terms. The electro- 
magnetic field is split in two fields, according to the dir- 
ection of propagation: one field for the photons in the 
forward direction (fc = 2), that of the stimulating laser 
and of the lost light, one field for the photons collected by 
the detector (fc = 1). In the rotating wave approximation 
we can take 



i?2 



\/7(l -P) 



The coefficient 7 > is the natural line-width of the 
atom, p is the fraction of the detected fluorescence light 
and 1 — p is the fraction of the lost light (0 < p < 1) 

pun mm]. 

b. Other dissipation terms. We introduce also the 
interaction with a thermal bath. 



i?3 — ^J^na^ , i?4 = \J^nGj^ , ?^ > 0, 
and a term responsible of dephasing (or decoherence) , 

k, > 0. 
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c. The laser wave. We consider a perfectly coherent 
monochromatic laser of frequency v > Q: 



fkit) 



in 



3fc2 



2\/7(l-p) 



(81) 



T is a time larger than any other time in the theory and 
the limit T — > +00 is taken in all the physical quantities. 
The quantity > is called Rabi frequency and Av = 
vq ~ v is called detuning. 

d. Master equation. With these choices the Li- 



ouville operator (151 becomes 



C{t)[p\ = - - [i/oaz + ilcr^t, p] + ^kd {(yzpo-z - p) 



+ 7(fT+l) (a_pa+--{P+,p} 



7n ( CT+po-- - - {P-,p] 



(82) 



The master equation ( 14 ) can be solved by using Bloch 



equations in the rotating frame [THl Sec. 8.2]. Indeed, we 

have 



T{t,s)[p] = e 



62 



(83a) 



C[p] = - - [v^a^ + nax, p] + "fkd {(Tzpaz - p) 



+ 7(n + l) (c7_p(T+--{P+,p} 



+ 7" ( CT+po-- - - {P-,p} 



(83b) 



The system reduced state turns out to be given by 
7y, = i{l + Kt)+iy(t)]e-V_ 



+ [x{t)-iy{t)]e~"''a+ + z{t)cjz] 



(84) 



where 



-At 



f (0) - 7 ■ 



^(i +71 + 2fcd) 

A = I 





A 



Av 
\-n + 2fcd) n 
-fl 7(1 + 2n). 



long time, the stimulating wave / and the local oscillator 
wave h must be produced by the same physical source 
and this means to take h proportional to /. So, by in- 
cluding any phase shift in the pase 1} already present in 
the definition (21), to describe homodyning we take 



hit) 



m)\ ■ 



With the choice (21 1 for /, we get h{t) 
The limit T - 



[Q,T 



(85) 



it). 



foo can be taken in Eqs. (63) and it is 
independent of the atomic initial state [44 . The result is 
\M Sec. 9.2.11 



T-y+00 
= 2TT'yp\s{d) •Xoq|^5(p), 



(86) 



oincl 
^hom 



where 



hm ^r'(M;^) 

T— i.+oo 



1 + 2p7 s{-d) 



A 



A^ + ^2 



A^) , (87) 




1) cosi? — , 



2^cq2/cq Siu 

sin i) — XcqVcq COS i9 




-jA~ 



Examples of inelastic homodyne spectra are plotted in 
Figure [1] for 7 = 1, n = fed = 0, p = 4/5. The Rabi fre- 
quency ri, the detuning Aiy and the phase z9 are chosen 
in order to get the deepest minimum of S^^^^^ in p = 2. 
Thus, in this case the analysis of the homodyne spec- 
trum reveals the squeezing of the detected light. Also the 
complementary spectrum is shown, in order to illustrate 
the role of the Heisenberg-type uncertainty relation ( 79 ) . 



One could also compare the homodyne spectrum with 
and without n and fed, thus verifying that the squeezing 
is very sensitive to any small perturbation. 



B. Heterodyning 



A. Homodyning 

The squeezing in the fluorescence light is revealed by 
homodyne detection, which needs to maintain phase co- 
herence between the laser stimulating the atom and the 
laser in the detection apparatus which determines the ob- 
servables Q(i;'i?, /i). To maintain phase coherence for a 



When the local oscillator and the stimulating wave are 
not produced by the same source, the phase difference 
cannot be maintained for a long time; in this case we 
have heterodyne detection. In the case of perfectly mono- 
chromatic waves, with a stimulating laser represented by 



(81 ), we take as local oscillator 



hit) 



l[0,T](i), I'lo^l^- 
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Figure 1. Sj,"''i(At; i?) with A;/ = 1.4937, ^ = 1.4360 and 
■d = -0.1748 (solid line), -d ^ \ - 0.1748 (dashed line). 



Again the limit T +oo can be taken in Eqs. (63) 
and it turns out to be independent of the atomic initial 
state and of i?. Let us set 



then, the final result is: 5*1101 (mS^Io) = 'S'hot(A'i '^lo) 

'het 



1 — + 



= + - v) + 5{^i + v)) 



= |K.„(m-^;0) + 5£L(m-^;V2) 



(89a) 



5hct (m; ^io) = ^lim 5r'(M; ^) - lpD{^i, v) 

1 — >-+oo 



+ 4 [SZnM - ^; 0) + SZliti - V- 7r/2) 

::'incl / ,, I ,. n\ i oincl 



+ SZL{^^ +v-Q) + SZL{^^ + ^; 7r/2)] , (89b) 



+ (/i + v)^ + (/i - v)^ 



t{0) 



A'^ + {^1+ w)2 A'^ + {^1- w)2 



(89c) 



Recall that s(0) = (1, 0, 0) and s{n/2) = (0, 1, 0). 

The inelastic heterodyne spectrum (89c) can also be 
written as 



S'hot ifJ-; V\o) = 1 + 2tTP [i;inol(M + V) + i;incl(M - v)] 



(90) 



-'incl 



(m) = ^ Rc (^(l,i,0) • [f(0) -if(7r/2)]) . 

(91) 



Properties of the heterodyne spectrum 

By explicit computations, it is possible to prove [161 
Proposition 9.3 and Remark 9.4 in Sec. 9.1.2] that 



Aiy = 



D{fi,v) = 



and that 



n = and fc^ = 



Di^i, v) = 0. 



In these cases the heterodyne spectrum reduces to a lin- 
ear combination of different homodyne contributions. 

a. The lower bound of the heterodyne spectrum. Be- 
ing SZt{p-,i'io) independent of any one of the two 
bounds in Theorem [l] implies 



(92) 



This means that it is impossible to see squeezing by het- 
erodyning. This is true not only in the model of this 
section, but in any physical set up for which the depend- 
ence on 'd is lost. 



2. The power spectrum 

Let us consider now the heterodyne spectrum as a func- 
tion of the frequency of the local oscillator in the case 
/i = 0. By particularizing the expressions (89), we get 



Qinc 
"-"het 



(0; ^^lo) = 2 (^;"om(^; 0) + SZ^{v; f )) + ^pDiO; v) 



= 1 



-t{Q)\-s{Q)- 



(93) 



^2 



„2 



It can be shown that S'hct(0; v\o) is proportional to the 
mean of the power of the heterodyne current |16l Sec. 
9.1.1]; so, as a function of t'lo, it represents the power 
spectrum. This interpretation can be strengthened by 
expressing iS'het(0; fc'io) in terms of the fields. Let us con- 
sider now the mode operators (66) in the case h{t) = 1: 

i-T 



1 



e''"*dBi(t) arip). 



These operators, together with their adjoint, satisfy the 
bosonic commutation relations ( 68 ) . By taking into ac- 



count that the iJ-dependent terms vanish in the limit 
T — > oo, from ( 76 ) we get 



<S'het(0; I'lo) = 1 + 2 lim Trr j n/(T)a^(z^io)aT(t'io) [ • 

T— f+oo I. J 

(94) 

So, the mean observed power spectrum is composed by 
the flat spectrum of the shot noise, conventionally set 
equal to 1, plus a term proportional to the mean number 
of photons in the temporal mode of frequency approxim- 
ately equal to i^io. 
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The fluorescence spectrum. The quantity 



5'het(0;i'io) - 1 7(= 



cq 



2/eq) 



Siv) 



-'incl 



iv) 



is interpreted as the fluorescence spectrum of the atom 
[TOl Sec. 9.1.2]; the normahzation is chosen in order to 
have its integral over v equal to the rate of emission of 
photons in the equilibrium state. For n = fc^; = this 
quantity coincides with the original MoUow spectrum [1^1 
Sec. 9.1.2.2], [m pp. 178-181], [HI P- 288]. In Figure [2] 
we give the plot of the inelastic part of the fluorescence 
spectrum in an asymmetric case {n ^ ^ 0) in which 
the Mollow triplet is well visible (fi large). 




Figure 2. 
fed = 0.7. 



Sinci(M) with Az^ = 1.5, Q, = 8.0, n = 0.01, 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have presented two formulations of 
the quantum theory of measurements in continuous time. 
The first one is based on the Hudson-Parthasarathy equa- 
tion, on quantum stochastic calculus and on the obser- 
vation of commuting field observables. The other one 
is based on the SSE, the stochastic master equation and 



the notion of conditional states. We have also shown how 
to connect the two formulations. Then, we have studied 
the properties of the observed output and in particular 
the output spectrum and the uncertainty relations on the 
spectra of incompatible quadratures. 

Being the two formulations of quantum continuous 
measurements equivalent, the output spectra can be de- 
duced from the classical SSE or from the quantum one, 
and the bounds (of Heisenberg type) on the spectra hold 
independently of the formulation. This point is conceptu- 
ally relevant, because this equivalence and the presence of 
uncertainty relations in both formulations show that the 
"classical" SSE is not a semiclassical approximation to 
some "quantum" theory, but it is itself fully "quantum" . 
However, let us stress that the proof of the bounds is 
based on the quantum field formulation and on the un- 
certainty relations for incompatible quadratures of the 
fields. On the other side we have shown that the most 
complete results on the probabilistic structure of the out- 
put can be obtained only in the stochastic formulation. 

The final part of the paper is devoted to a relat- 
ively simple example which allows to show squeezing in 
the homodyne spectrum and to illustrate the differences 
between homodyning and heterodyning. 

Various generalizations of the theory presented in this 
paper are possible, first of all by introducing direct de- 
tection ^iliZlilllSlElHaaiSO] and Markovian feedback 
[ini HB [ini SI m IMSS] . The quantum trajectory ap- 
proach can be generalized to introduce also feedback con- 
trol with delay, coloured noises, various non-Markovian 
efi^ects [50H55] . The simplest non-Markovian contribu- 
tion is to take the laser wave / to be random [531 123 . 
This means to consider as initial state of the field a clas- 
sical mixture of coherent states, the mean of grif)- Even 
the local oscillator wave h can be random, which again 
means to take a mixture of coherent vectors as state of 
the local oscillator [551 ■ This randomness can be used 
to introduce more realistic models of laser light, not only 
perfectly coherent monochromatic waves, but also waves 
exhibiting some coherence time. This allows for a bet- 
ter analysis of the differences between homodyning and 
heterodyning |55j . 
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